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"^s ' Abstract 

Transaction costs appear in financial markets in more than one form. There are several results in the 

literature on small proportional transaction cost and not that many on fixed transaction cost. In the 

present work, we heuristically study the effect of both types of transaction cost by focusing on a portfolio 

optimization. Here we assume the presence of fixed transaction cost and that there is a balance between 

fixed and proportional transaction cost, such that none of them dominates the other, asymptotically. We 

find out that the deviation of value function, when the fixed transaction cost is e, from the Merton value 
• , i 

function, without transaction cost, is of order £2 which is different from the pure proportional cost of 

£3 . Based on this, we propose an expansion for the value function in terms of powers of £2 . 

1 Introduction. 

Trading strategies developed by Atkinson, Pliska, and Wilmott [I] consider the effects of a small fixed cost 
payed by the investor independently of the volume of the transaction. The case of small cost proportional 
to the volume of the transaction has been studied by Shreve and Soner [7J and Goodman and Ostrov [3J. 
In both cases the emphasis is on the asymptotic behavior of the trading strategies as the transaction cost 
vanishes and the problem approaches the idealized no transaction cost problem studied by Merton [3]. In 
this paper we consider the case when the investor faces both fixed and proportional vanishing transaction 
costs. We find the balance between the two costs that makes no cost dominant and describe the asymptotic 
optimal trading strategy as a function of the proportional cost parameter, A, and the fixed cost parameter, 
e. 

The case of small proportional cost has been extensively studied. Shreve and Soner in [7J and Janecek and 
Shreve in [3] present rigorous arguments based on viscosity solutions to show that in the case of power utility 
and a single stock, the value function of the investment-consumption problem under small proportional cost 
A converges to the Merton value function as fast as 0(Aa ). Rogers [6 1, proves the same result using a more 
probabilistic argument. For general utilities and multiple stocks, in [5], the authors provide the result using a 
heuristic argument based on a perturbation method. Goodman and Ostrov [5J show the same result based on 
a dual probabilistic argument which also provides a clear understanding of the quasi-steady state probability 
density of the optimal portfolio. The traditional asymptotic expansion for the value function , /, is given by 
the equation 

/ = f° + A 1 / 2 + A/ 3 + \\f 4 + o(A*) • (1.1) 

In the recent paper [5], the authors give a rigorous proof of the equation \f — f°\ ~ 0(A3) for general 
utilities under some assumptions on the regularity of Merton value function and on the Merton ratio. Their 
argument are based on the theory of viscosity solutions. 
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There are fewer studies in the case of small fixed cost. Atkinson, Pliska and Wilmott [T] consider the 
effect on trading strategies of a small fixed cost payed by the investor. They perform a formal asymptotic 
analysis of the value function associated with the problem at hand and obtain an expansion of the form 

/ = /o + £/ 4 + o(e) . 

We study the case of a single stock. The argument follows the heuristics used by Goodman and Ostrov [2] 
where the effect of the costs is decomposed in trading cost and oportunity cost, which is due to deviations from 
the idealized portfolio. The (approximate) optimal strategy keeps the portfolio position inside a hold region, 
W, centered around the idealized portfolio position. A transaction takes place only when the portfolio's 
position is on the boundary of T-L to move the portfolio position to the boundary of a destination region, 
T>, also centered around the idealized portfolio position and contained in H. Assuming that the destination 
region is fixed, an expansion of the hold region decreases the trading cost but increases the opportunity cost. 
On the other hand, if the hold region is fixed, an expansion of the destination region will certainly increase 
the opportunity cost but it is not a priori clear how it affects the trading cost. The optimal trading strategy 
is obtained by quantifying and balancing these effects. The main tool the is equilibrium probability density 
for the portfolio position, u, which turns out to be constant in T> and linear in T-L \ T). 

The heuristic analysis has been shown to be the dual equivalent of the perturbation theory applied to 
/, the value function for the expected utility of the portfolio at a final time T, in the case of proportional 
transaction cost. In section [2] wc perform both the heuristic analysis and the perturbation analysis of the 
value function / in the fixed transaction cost setting as a sanity check. In section [3] we rely only on the dual 
heuristics when both fixed and proportional costs are present. When the hold (destination) region has a 
radius equal to 7 (77) the opportunity cost is proportional to 7 2 + n 2 and the transaction cost is proportional 
to 2 _ 2 e ■ Optimizing with respect to 7 and r\ shows that the balance is obtained when e 3 — A 4 and in 

this case 7 is of order e 1 / 4 and 7 — 77 is of order e 3 / 4 . The balance between the transaction cost parameters 
A and e is the only one such that both costs have the same effect on the utility function /. 

2 Small fixed transaction cost. 

2.1 Idealized no transaction cost problem 

The classical allocation problem for an investor with z dollars invested in a portfolio at time t is to maximize 
the expected utility U of the portfolio dollar worth at a final time T; that is to find 

f(z,t) = supE Ztt [U(Z(T))}. (2.1) 

Here the supremum is taken over all the admissible strategies available to the investor by buying or selling 
a risky asset. This problem was studied by Merton in [4]. We recollect the main results in order to fix the 
notation, which is borrowed from [2]. The investor's portfolio is composed of X dollars in a stock whose 
price in dollars S follows a geometric brownian motion 

dS = )j,Sdt + aSdB, 

and a money market account worth Y dollars and growing at a constant risk free rate r. The SDE's for the 
stock and money market dollar worth are 

dX = fiXdt + aXdB (2.2) 

dY = rYdt, (2.3) 

respectively and the SDE for the portfolio dollar worth Z = X + Y is 

dZ = [(/x - r)X + rZ]dt + aXdB . (2.4) 



We can write down the SDE solved by / = f(Z,t) as 

df = I f t + [(ji - r)X + rZ]f z + ^a 2 X 2 f zz \dt + oXf z dB . (2.5) 

The investor has the freedom to transfer money at anytime between the stock portion of his portfolio and 
the money market account without incurring in a transaction cost. The optimal strategy for the investor 
will depend only on the the time t and the portfolio dollar worth Z; this is X = x(Z,i). The Hamilton- 
Jacobi-Bellman equation for the optimization (|2.1j) is 

= supj/t + [(n - r)x + rZ]f z + ~a 2 x 2 f„\ . (2.6) 



The solution is the Merton ratio 



x = m{Zit) = Jj^M±_ (2.7) 

°~ Jz 



I zz 



Under the optimal control X we have f(z,t) — E Z ^\U(Z(T))] = E z j[f(Z(T),T)] and therefore the process 
f(Zt,t) is a martingale, which we write in differential notation as 

E[df] = . (2.8) 

Under a suboptimal control X we will obviously have E[df] < 0. The martingale property (|2.8[) together 
with equation (I2.5[) show that the value function / solves the PDE 



l (M-?f/ g 2 

2 <J 2 f zz ' '"" " "" (2.9) 

f{z,T) = U{z) 



ft ~ o 27 + rzf* =0, for t < T 



In the rest of this paper we will denote the solution of the previous nonlinear equation by f° and make clear 
that it is the value function under zero transaction cost. 

2.2 Heuristic analysis of small fixed transaction cost. 

We will assume that the investor's broker charges a small fixed cost of e dollars for buying or selling any 
number of stock shares at a fixed time. This model is only accurate for very liquid stocks and small changes 
in the stock portfolio. In this case the flat fee is the main source of transaction cost. Borrowing notation 
from [2], let L(t) be the dollar amount of cash spent buying stock up to time t, and M(t) be the dollar worth 
of all stock sold up to time t. These are the controls available to the investor. The SDE's for the stock and 
cash components of the investor portfolio become 

dX = Lt,Xdt + <jXdB + dL-sS dL>a dt-dM (2.10) 

dY = rYdt + dM - eS dM>0 dt - dL . (2.11) 

Here 5dL>o is the sum of delta functions at the times where dL > and an analogous definition holds for 
SdM>o- At this point it is standard to perform the change of variables £ = X — m(Z,t) and work with £ 
instead of X. Here m(Z, t) is the Merton ratio defined by equation (|2.7[l . This is natural because the optimal 
X in the case with no transaction cost is X — m(Z, t). Using Ito's formula, the SDE's for Z = X + Y and 
£ = X — m are 



and 



dZ = [(p-r){m + £)+rZ}dt + o-{t + m)dB-e(5 dL> odt + 5 dM >odt) (2.12) 



d£, = \jj,(£ + m) -m t - -m zz <j 2 (f + to) 2 - m z (((i - r)(m + g) + rZ)]dt 



+ (1 - m z )o-(£ + m)dB + A£ (5 dM >o + $dM>o)dt 



The symbol A£ denotes the jump in the value of £ at time t. The allocation problem becomes 

f(z^,t)=snpE ZiLt [U(Z(T))}. 

Our task is to perform an asymptotic analysis for the value function / around the zero transaction cost 
value function / defined in subsection 12.11 In what follows we will give heuristic arguments that will be 
justified later in subsection 12.31 We start by considering strategies such that there is a small region of 
the form £ 6 (—j(Z,t),j(Z,t)) where no transaction takes place and when £ touches the boundary =F7 we 
perform a transaction in order to make £ = 0. This makes sense because the transaction cost is independent 
of the transaction amount, so the optimal strategy must make £ ~ immediately after £ = ±7. 

The study of the dynamics of £ will be the main tool in the asymptotic analysis. The process £ has 
behaves like a diffusion before hitting the boundary ±7 and it is restricted to a confined region. All the 
possible values of £ are visited in a small time interval and therefore £ must reach statistical equilibrium 
much faster than Z. In other words, in all subsequent calculations we will think of Z and t as constants and 
£ will be assumed to be in equilibrium with density «(£). The leading order of the SDE (|2.13j) satisfied by 
£ is calculated by first dropping the terms that contain an e factor. We can also drop terms that have a £ 
factor because |£| < 7 and we are assuming that 7 is small. Finally, we drop the drift terms because the 
process stays inside a small region. Thus, up to leading order, 

d£~adB + A£, (2.14) 

with a = (j(l — m z )(£ + m). Under such an strategy, the differential equation for the equilibrium probability 
density w(£) is 

1 a 2 

-a 2 Uii + —S = (2.15) 

ut = T^r at £ = ±7. (2.16) 

7 

The constant a 2 /^ 2 is the rate of particles hitting the boundary £ = ±7 per unit time. When they hit the 

2 
boundary, they are transported to £ = and hence we obtain a source term proportional to ^ at £ = 0. 

The solution is 



with Newmann boundary conditions 
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The optimal boundary 7 is obtained by maximizing E[df(Z, £, t)] over 7. Instead, we will maximize 
E[df°(Z,t)]. This approximation will be justified later in subsection 12. 31 . Ito's formula and equations (J2.12J) 
and (|2.9j) show that the idealized Merton value function process f°(Z,t) solves the SDE 



(p ~ r)!if Q z + ^mif zz + ^ff zz + Af(5 dL>0 + S dM>0 ) 



dt 



df = 

+ a{m + 0f°dB 
We will calculate the £- equilibrium expected value of the previous equation. Use that in equilibrium 

E[£] = 



and 



j — 1 



-1 
1 2 

6 7 



to get 



E[df] 



2 2 

12 



a 2 f° zz + Af(E[5 dL>0 ]+E{5, 



dM>0\ 



(It . 



When dL > or dM > the process Z jumps from Z to Z — e and therefore A/ ~ — ef z . Thus, 



£[d/°] 



12 



/«-4"(£feo]+£[Wo]) 



dt. 



(2.17) 



In order to calculate E[5dL>o], apply Ito's formula to the process £ 2 to obtain the equation 

d(£ 2 ) =s 2£adB + a 2 di + A(f)S dL>0 dt + A(£ 2 )<5 dM >o^ . 

The expected value of the left hand side of the previous equation is approximately zero because £ is assumed 
to be close to statistical equilibrium, thus 

~ a 2 dt - ~f 2 E[6 dL>0 }dt - ~f 2 E[6 dM >o}dt . 

The symmetry of the boundary of (—7, 7) gives 

E[5 d L>o] 



1 



2 7 2 ' 
and plugin in the last equation in (|2.17|) we obtain an approximation for 



E[df°] 



dt 



in terms of e and 7: 



E[df°] 
dt 



2 2 ? 
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Maximizing the right hand side of the previous equation over 7 gives 

= ! -12- 



n 2f0 \ 1/4 

a Jz \ 1/4 



T 2f0 

/ zz 



and therefore 



E[df°] 



o 2 l. 



.1/2 



(2.18) 



dt \ 12 

We conclude that the fixed cost e shifts the value function by e 1 / 2 . 

2.3 Fixed cost asymptotic analysis. 

The heuristic arguments of the previous section suggest that £ moves in a region of size e 1 ' 4 . We define a 
rescaled variable £ = e _1 / 4 £ but we drop the tilde to avoid cumbersome notation. We obtain from equation 
(|2.13p that the rescaled process £ solves the SDE 



di =e- l/i [n(e l/i £ + m) - m t - lm zz <j 2 {e 1/4 £, + to) 2 - m z ((p - r)(m + e 1/4 + rZ)]dt 

+£- 1 / 4 (l - m z )<r{e- 1/4 £ + m)dB 
+AC- 



(2.19) 



Our control consists of a continuation region of the form (— j3, 7) and two points rj, 6. When the process £ 
hits the boundary £ = 7 the investor buys stock in order to get £ = r\. Similarly, when the process £ hits 



the boundary £ = — /3 the investor sells stock in order to get £ = —8. In the continuation region we have 
dL = dM = 0. Using this fact, the SDE's (|2.19p and (|2.12[) and Ito's lemma we conclude that the process / 
solves Bellman's equation (in the continuation region): 



0=£- 



1/4 He 1/4 £ + m)-m t - Im zz a 2 (e 1/4 £ + ™) 2 - m z (( M - r)(m + e 1 '^) + rZ)]fr 



2 



+ l -e-^{l-m z fa\^^ + m ff a 



2 

+ £ - 1 / 4 (l-m z )a 2 ( £ 1 / 4 e + m) 2 ^ (2 ' 20) 

+[(/* - r)(m + e 1 ^) + «]/, + I<r 2 (e 1/4 £ + m) 2 f zz 

+ff 

On the boundary of the hold region the jump terms in equations (|2.19|) and (|2.12[) dominate the Hamilton- 
Jacobi-Bellman equation so we must have 

f(€,z,t)=f(r),z-e,t), C = 7 (22l) 

f(£,z,t)=f(-e,z-E,t), £ = -/3 

The optimality boundary conditions (smooth pasting, take <9 7 and d n of the conditions at the boundary) are 

m,z,t)=o, e = 7,^ (222) 

We propose an asymptotic expansion in powers of e 1 / 4 because this is the power of e that we used to 
rescale £. The asymptotic expansion reads 

/ = /»(*, t) + eWfHS, z, t) + e 1 / 2 / 2 ^, z, t) + e 3 / 4 / 3 (£, z, t) + e/ 4 (£, z, t) + o(e) (2.23) 

The function f° is the value function in the e = case. Using the boundary condition (|2.21|l and performing 
a Taylor expansion around (£, z) we obtain 

0=-ef z + O(e 2 ) + fs(r 1 -£) + •■• . 

This means that f z is 0(e) smaller than /j and therefore J 1 , / 2 and / 3 do not depend on £. Now we use the 
expansion (|2.23|) and collect e terms in the boundary conditions (|2.21j) and obtain 

f 4 (l,z,t) = -f z +f 4 ( V ,z,t) 
f(-0, z ,t) = -f z +f\-8,z,t) 

Collecting e terms in the optimality boundary conditions gives 

fl(£,z,t) = 0, £ = 7,-/3,r/,-0. (2.25) 

Finally, we plug in the asymptotic expansion in Bellman's equation (|2.20[) . The 0(e) equation involves only 
/° and it is just the Merton equation 

= /° + [(/, - r)m + rz]f z + \a 2 m 2 f zz . 

The (^(e 1 / 4 ) equation is 

0=[(p-rnf z +a 2 ^mf° zz 
+ f} + [{jjl - r)m + rz]fl + a 2 m 2 f l zz . 



Since to is equal to the Merton ratio the terms involving / in the previous equation cancel each other 
and we conclude that f 1 solves Merton's equation. The final condition is f 1 (z,T) = because f°(z,T) = 
f(z,£,T) = U(z) and therefore, by uniqueness, f 1 (z,t) = 0. The 0(e 1 ^ 2 ) equation is 

= f 2 + [( M - r)m + rz].f z + W f zz + \a 2 i 2 f zz + 5(1 - m,)Vm 2 /| s 

z z L (2.26) 

= K + Bi 2 + -a 2 f^ , 



where 



2 



K(z, t) = f* + [(/i - r)m + rz]f 2 + -a 2 m 2 f 2 z 

B{z,t) = \a 2 f zz 
a{z, t) = (1 — m z )am . 
Integrating with respect to £ twice we obtain 

where D(z,t) and C(z,t) are constants of integration. Now we have the six equations given by (|2.24|) and 
(|2.25j) and the six unknowns 7, ft, ??, 9, C and D. The unknown D does not appear in the equations and can 
be set to any value that we like. Modulo D, the unique solution is C = r\ = 9 = and 

f a 2 f°\ 1/4 

3 Fixed and proportional transaction cost 

Assume that the fixed transaction cost is $e for buying or selling any number of stocks and the proportional 
cost is X% of the transaction in dollars. We will consider strategies such that £ = ±77 immediately after 
£ = ±7. The new SDE's solved by X and Y are a small modification of equations (|2.10|) and (|2.11j) . namely 

dX = f iXdt + aXdB + (l- X)dL - sS dL>0 dt - dM (3.1) 

dY = rYdt + (1 - X)dM - eS dM >odt - dL . (3.2) 

Using Ito's formula, the SDE's for Z = X + Y and £ = X — m are 

dZ = [(/! - r)(m + 0+ rZ]dt + <r(£ + m)dB - X(dL + dM) - e{S dL>0 dt + S dM>0 dt) (3.3) 

and 

d£ = [m(£ + m)- m t - ^m zz a 2 {£, + to) 2 - to 2 ((/i - r)(m + £) + rZ)]dt 

+ (1 - m z )cr(£ + m)dB + A£ 

Applying the same heuristics used in section [2] we will think of Z and t as constants and £ will be assumed 
to be in equilibrium with density u(£). The leading order of the SDE (|2.13[) satisfied by £ is calculated by 
first dropping the terms that contain an e or A factor. We can also drop terms that have a £ factor because 
|£| < 7 and we are assuming that 7 is small. Finally, we drop the drift terms because the process stays inside 
a small region. Thus, up to leading order, 

d£ ~ adB + A£ , (3.5) 



with a = cr(l — m z )(^ + m). Under such an strategy, the differential equation for the equilibrium probability 
density u(£) is 



1 



fl2u « + Tu~2 Z2\ S V + 77772 -5T<S-, 



2" -« ' 2( 7 2 -77 2 ) " 2( 7 2 -77 2 ) 
with Newmann boundary conditions 

1 







u i = T- 2 



at e = ±7 ■ 



(3.6) 



(3.7) 



The constant £_ a is the rate of particles hitting the boundary £ = ±7 per unit time. This fact can be 
shown, as in section [2 by applying Ito's formula to the process £ 2 , namely 

d(£ 2 ) ~ 2iadB + a 2 dt + A(£ 2 )S dL>0 dt + A($ 2 )6 dM >odt . 



Taking the equilibrium expected value of the previous equation and using the symmetry of the boundary 
we obtain E[SdL>o] = 2 ( °- 2 ) • When the particles hit the boundary, they are transported to £ = ±77 and 

at £ = ±7;. The solution is 



2( 7 2-rP) 

hence we obtain a source term proportional to 



«(0 



1 nt\-v 



7 + V V7 2 _ ff 



The optimal boundary 7 and location 77 are obtained by maximizing E[df(Z, £,£)] over 7 and 77. Instead, 
we will maximize E[df°(Z, t)]. Ito's formula and equations (12.12)) and (|2.9p show that the idealized Merton 
value function process f°(Z,t) solves the SDE 



df° = 



0* - r)£/° + v 2 mtf° z + ^-f° zz + A/°(<5 dL>0 + <Wo) 



dt 



(3.8) 



+ °{m + Of z dB 
We will calculate the £- equilibrium expected value of the previous equation. Use that in equilibrium 

E[£] = , (3.9) 

and 

m 2 } = J" e a «(o# = \(h 2 + v 2 ) , (3.10) 

to get 

V 2 (7 2 + ?7 2 ) 



E[df] 



12 



2 f° zz + Af (E[6 dL>0 ]+E[6, 



dM>0\ 



dt. 



When dL > or dM > the process Z jumps from Z to Z — e — X\dX\. Since X = £ — r7i(Z, i) when a jump 
takes place we have \dX\ ~ |d£| =7 — 77 and therefore A/ ~ — [e + A( 7 — T])]/®. Using that both E[5 d L > o] 
and -E[cW/>o] ar e equal to 2 , °_ ^ we obtain the equation 



£ ' /r •■C(7 2 +^) + A((T-l)A + E ), 



f/f 



^,2 _ ? ^2 



(3.H) 



er 2 f° 

where C = 12 " and it = — /. 



z ° 2 °_ 2 . We use the change of variables m = 7 + 77, n = 7 — 77 and write 






F(m,n) = —(m 2 +n 2 ) H (nA + e) 

2 777,77 



In order to find the optimal values of m and n we calculate 

F m = Cm — (nX + e) 

m z n 

F n — Cn K e 



and 



D'F = ( ^ ' c Y^R 



Using that the gradient of F is zero at the optimal controls we obtain 

C 2 n 8 {nX + e) = R 2 e 3 
Now write n = cue/X in order to get 



C 2 u\lo + 1) = R 2 
e 3 = A 4 



R rl/4 (3-12) 



Clo 3 



£ 



n = cV/ 4 



We undo the change of variables and conclude that 

R £ l/4 +w 3 £ 3/4 

■y — S=diL 

' 2 

ry = 7 - w 3 e 3/4 . 

Notice that e 3 — A 4 is the only balance between the two types of transaction cost that makes the shift of the 
value function under proportional transaction cost equal to e 1 ' 2 , which is the shift due to the proportional 
transaction cost. 
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